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The dynamics of the kicked-rotor, that is a paradigm for a mixed system, where the motion in 
some parts of phase space is chaotic and in other parts is regular is studied statistically. The evolution 
( Frobenius-Perron ) operator of phase space densities in the chaotic component is calculated in 
presence of noise, and the limit of vanishing noise is taken is taken in the end of calculation. 
The relaxation rates ( related to the Ruelle resonances ) to the invariant equilibrium density are 
calculated analytically within an approximation that improves with increasing stochasticity. The 
results are tested numerically. The global picture of relaxation to the equilibrium density in the 

^^ ' chaotic component when the system is bounded and of diffusive behavior when it is unbounded is 

^\ I presented. 
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Q ! Statistical analysis is most appropriate for the exploration of global properties of systems that 

. exhibit complicated dynamics l|l]-ffl . For chaotic systems, ?7, the evolution operator of distributions 
00 • 

CN ' of phase space trajectories that is sometimes called the Frobenius-Perron (FP) operator describes 

T-H . the statistical properties of the dynamics. For many idealized systems exponential relaxation to the 
> ■ 

^\ ' equilibrium density takes place. This was established rigorously for hyperbolic systems (A systems), 

m : __ 

O ■ like the baker map [0-|§]. The relaxation rates related to the Ruelle resonances that are poles of the 

O ' . . 

T-H ' matrix elements of the resolvent R = {z — U)^-^ in a space of functions that are sufficiently smooth 

0\'. 

0\ . 1^. These poles are inside the unit circle in the complex z plane while the spectrum of U is confined 

H ' to the unit circle because of unitarity. Most physically realistic models are not hyperbolic, but mixed 

Y ■ where the phase space consists of chaotic and regular components. For mixed systems sticking to 

c^ ' regular regions takes place. If the regular regions are small this effect is negligible for finite time, 

Q that may be long, much longer than the time relevant to the experiment. In this letter the FP 



operator and the relevant approximate relaxation rates are calculated analytically and numerically 
for the chaotic component of the kicked rotor, that is a mixed system 0. 

The kicked rotor is a paradigm |^ for chaotic behavior of systems where one variable may be 
either bounded or unbounded in phase space. If it is unbounded diffusion is found for the classical 
system PJ^ . In quantum mechanics this diffusion is suppressed by a mechanism similar to Anderson 
localization |1^. The kicked rotor is defined by the Hamiltonian 

n = ^J^ + KcoseY,Sit-n), (1) 

^ n 

where J is the angular momentum, 6 is the conjugate angle (0 < ^ < 27r) and K is the stochasticity 
parameter. Its equations of motion reduce to the standard map 6 = 6 + J and J = J — Ksin6, 
where {6, J) and {6, J) are the angle and the angular momentum before a kick, and just before the 
next kick respectively. For K > K^^ 0.9716 diffusion in phase space was found. 



In the present paper the FP operator will be calculated for the kicked rotor on the torus: 
(0 < J < 27TS ;0 < 6 < 2tt), where s is integer. The operator is defined in the space spanned by the 
Fourier basis: 

11 / kj\ 

(pkm = {J6\km) = —r=—r== exp{ime) exp i — . (2) 
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The FP operator was studied rigorously for the hyperbolic systems and many of its properties are 

known [0-0] • It is a unitary operator in C^, the Hilbert space of square integrable functions. 

Therefore its resolvent 

1 1 °° . 

R{z) = ^ = -Y: U-z- (3) 

z — U z ^^Q 

is singular on the unit circle in the complex z plane. The matrix elements of R are discontinuous there 
and one finds a jump between two Riemann sheets. The sum (^ is convergent for \z\ > 1, therefore 
it identifies the physical sheet, as the one connected with this region. The Ruelle resonances are the 
poles of the matrix elements of the resolvent, on the Riemann sheet, extrapolated from |z| > 1 [Q. 
These describe the decay of smooth probability distribution functions to the invariant density in a 
coarse grained form 0. In spite of the solid mathematical theory there are very few examples where 
the Ruelle resonances were calculated for specific systems P,|l. For the baker map it is easy to see 
that as the resonances approach the unit circle, corresponding to slower decay, they are associated 
with coarser resolution in phase space 0. For the kicked rotor (|l]) the operation of ?7 on a phase 
space density p is Up{0, J) = p{6 — J, J + Ksin{6 — J)). To make the calculation well defined noise 
is added to the system. If noise that conserves J and leads to diffusion in 6 is added to the free 
motion, the matrix elements of JJ in the basis (0) are: 

{k2m2\tj\kimi) = J^a-mi I — — j exp ( — —ml j 6k2-ki,m2s- (4) 

For a = the operator is unitary as required. It is shown explicitly that addition of the noise acts 
effectively as coarse graining and the resulting evolution operator is not unitary (see also ||11|| ). 
For large stochasticity parameter K, it is shown here that in the Fourier basis (0) the slowest 
relaxation modes, in the limit of infinitesimal noise, are found to be identical to the modes of the 
diffusion operator |^. Also the fast relaxation modes are calculated analytically in the present 
work, and the approximate analytical results are tested numerically |^. These modes are not 
related to the spectrum of the FP operator that is confined to the unit circle. We believe we found 
all relaxation rates for distribution functions that can be expanded in terms of the basis functions 
(I). The immediate question is how is it possible that this description, that was established only 



for hyperbolic systems, holds for a mixed system. It is clearly approximate, and holds for large 
values of the stochasticity parameter K, since then most of the phase space is covered by the chaotic 
component. The physical reason for the decay of correlations is, that in a chaotic system, the 
stretching and folding mechanisms lead to a persistent flow in the direction of functions with finer 
details, namely larger \k\ and |m| in our case. Consequently the projection on a given function, for 
example one of the basis functions (0) in our case, decays [|T^ . The crucial point is that this function 
should be sufficiently smooth. This argument should hold as an approximation also for the chaotic 
component of mixed systems. For smaller values of K the weight of the regular regions increases. In 
such a situation, in the limit of increasing resolution the resonances related to the regular component 
are expected to move to the unit circle in the complex z plane, corresponding to the quasi-periodic 
motion, while the resonances associated with the chaotic component stay inside the unit circle [|13|. 
How is the FP operator related to the quantum mechanical evolution operator? It was shown 
numerically for the baker map that if both operators are calculated with finite resolution they 



exhibit the same Ruelle resonances ||TT| . Noise and coarse graining are introduced in field theoretical 
treatment of chaotic systems [|l4| , pr5| . Since the FP operator plays an important role in these theories 
our work is of relevance there. It also justifies some assumptions made in the calculation of the typical 
localization length for the kicked rotor p!6| , p^ . 

We turn now to calculate the Ruelle resonances for the kicked rotor with the help of the evolution 
operator (^). The calculation will be done for finite noise a and then the limit a — > will be taken. 
These are the poles of matrix elements R12 = {kimi\R{z)\k2m2) of the resolvent operator R of (^ 
when analytically continued from outside to inside the unit circle in the complex plane. It is useful to 
introduce R'{z) = 1/(1— zU) = Z^^o-^"^" ^^^^ i^ convergent inside the unit circle, because \\zU\\ < 
1. Its matrix elements are R^2 = {kirriilR' {z)\k2m2) = Z^^o'^n-^"' where a„ = (A;imi|f/"|A;2"^2)- The 
relation between the matrix elements inside and outside of the unit circle implies that if Zc is a 
singularity of R12 then 1/zc is a singular point of R^2- Consequently the first singularity of the 
analytic continuation of Ri2{z) from inside to outside the unit circle gives the first singularity one 
incounters when analytically continuing Ruiz) from outside to inside the unit circle, i.e. it is just 
the leading nontrivial resonance. It is determined from the fact that it is the radius of convergence 



r of the series for R'^2 is given by the Cauchy-Hadamard theorem: r ^ = lim„^oo sup \/\an\ |jT8[ 



The calculation of the coefficients a„ = {k0\tj"'\k0) is performed using the resolution of the identity 
( introducing intermediate \kimi){kimi\ ), and then substitution of (Q) and summation over the /cj 
leading to: 



mi m2 m„_i ^=1 \ "^ / 

where vfiQ = vfin = while M,^ = I]'=o"^« ^^^ ^F — "^z-i ~ ""^i- The calculation is performed for 
large s and K and the limits are taken in order |]^: 

(1) s ^ oo, (2) K ^ oo, (3) a^O. (6) 



For a sufficiently low mode so that < kK/s « 1, the leading order term in K/s and l/yK is 



an ~ 



fc2ir2 



1-^(1- 2J.(ir)e-^^ 



4s 
The resonance closest to the unit circle, Zk is the inverse of the radius of convergence. Here Zk 

^-(k^/s^)D{K)^ with 

D{K) = ^{l-2UK)e-\ 



(7) 



that is just the value of the diffusion coefficient D{K) found in [|T9|. In the limit of a ^ these 
are the relaxation rates in the diffusion equation. The analysis of the off-diagonal matrix elements 
Qn = [kmlU^-lk'm') leads to the same result. 

In order to obtain the fast relaxation rates we have to calculate matrix elements that do not exhibit 
slow relaxation, because such relaxation if present dominates the long time behavior. For this purpose 
we calculated the relaxation rates of disturbances from invariant density that involve functions from 
the subspace { 10,771) } with m 7^ and calculate a„ = {Om\tJ'"'\k'm'). Again the resolution of the 
identity is introduced ( introducing intermediate \kimi){kimi\ ), and summation over the ki yields 
a non vanishing result only if k/s = q that is an integer. The expression for a„ is found to be 



independent of s. The resulting resonances ( for large K ) are Zp = Jj2p{pK) exp {—a'^p'^/2), where 
p = \m\ OT p = \q\, where q = q ii q ^ and q = m' ii q = 0, depending which choice gives the 
larger absolute value 0. For q = m' = 0, and ??i 7^ one finds a„ = 0. If rri = m' = g = the 
only contribution is when all rrii vanish and then a„ = 1 for all n, resulting in the resonance z = 1, 
corresponding to equilibrium. 

The FP operator is the evolution operator U in the limit of vanishing noise. Therefore the Ruelle 
resonances are the poles of matrix elements of the resolvent R in this limit. They form several 
groups. There is 2:0 = 1, that is related to the equilibrium state. The resonances corresponding to 
the relaxation modes related to the diffusion in the angular momentum are: 

Zk = exp[-^^{l-2UK))\. (9) 

The resonances related to fast relaxation in the 9 direction are: 



J2r>{pK). 



(10) 



In certain cases this result does not hold for small intervals around K* . The relaxation rates are 
7fc = \lnzk\ and % = |/n|2p||. 

The analytical results that were obtained as the leading terms in an expansion in powers of 
1/yK were tested numerically for finite K and cr = 0. For this purpose the correlation function 
^fgi''^) ~ {f\U"'\9) "^^s calculated numerically. For distributions g and / from the Fourier basis (0), 
projected on the chaotic component, the relaxation rates are expected to take the values 7^ or 7^. 
For the diffusive modes one expects 7^ = {k'^/s'^)D{K), where D{K) is the diffusion coefficient (|]) 
with a = 0. The values of D{K) were extracted from this relation for various values of k and s and 
presented in Fig. 1. 
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FIG. 1. The difFusion coefficient D as extracted from the relaxation rates: the first mode (squares) , the second mode 
(stars), the fifth mode (circles) and off diagonal correlation functions (triangles), compared to the theoretical value (solid line). 
The dashed line represents the approximate error. The values of D obtained by direct simulation of propagation of trajectories 
are marked by diamonds. 



For large values of K, excellent agreement with the theory is found: the value of D is found to 



be independent of k and s and it agrees with (§). For relatively smaller values of K, the value of 
diffusion coefficient for some values of K is larger than the one that is theoretically predicted. The 
theoretical errors were estimated from the next term of the formula of Rechester and White for the 
diffusion coefficient |T^. In order to observe the rapidly relaxing modes the correlation function 
Cfg was calculated for g = 0^^,' so that q = k/s is an integer and / = 0om- In Fig- 2 the numerical 
estimate for 7^ is compared with the theoretical prediction obtained from (|TU|). 
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FIG. 2. The fast relaxation rates 7^ for / — (j>o\,g — <j)02 (triangles), compared to the theoretical value (solid line). The 
dashed lines denote the theoretically estimated error. Here we used s = 1 and A'^ — 10*. 



The error in the theoretical prediction is estimated as the value of the next order contribution to 
a„. The main reason for disagreement between the theory and the numerical simulations is sticking 



to the islands of stability and accelerator modes ^0 



Finite noise leads to the effective truncation of the evolution operator (^ . In the basis (^ it 
means that it results in limited resolution. Moreover for a > the operator U is non unitary. The 
approximate eigenvalues of If given by @ that were found in this work are 1, z^ and Zp. Because 
of the effective truncation, tp^, the eigenfunction of U, can be expanded in terms of the basis states 
(0). The relaxation rates of these eigenstates are —ln{zk) and —ln{\zp\). In the limit a ^ the 
evolution operator is unitary, and ip^ approach some generalized functions while z^ and Zp approach 
the values (p|,p!0[). These are the Ruelle resonances similar to the ones found for hyperbolic systems 
such as the baker map 0. Here noise was used in order to make the analytical calculations possible. 
In real experiments some level of noise is present, therefore the results in presence of noise are of 
experimental relevance. 

In summary, the Ruelle resonances, that were found rigorously for hyperbolic systems can be 
used for an approximate description of relaxation and transport in the chaotic component of mixed 
systems. The relaxation of distributions in phase space to the invariant density takes place in stages. 



First the inhoinogeneity in 9 decays with the rapid relaxation rates 7p and then relaxation of the 
inhomogeneities in the J direction takes place with the relaxation rates of the diffusion equation. In 
the limit s ^ oo the inhomogeneity in 9 relaxes and then diffusion in the momentum direction takes 
place. 
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